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Abstract 

We employ the variational method for the Sturm-Liouville eigenvalue problem to analytically 
investigate the properties of the holographic superconductors. We find that the analytic method is 
still powerful when the backreaction is turned on. Reducing step size in the iterative procedure, 
we observe that the consistency of results between the analytic and numerical computations can be 
further improved. The obtained analytic result can be used to back up the numerical computations 
in the holographic superconductor in the fully backreacted spacetime. 
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I. INTRODUCTION 



As a powerful tool to analyse strongly coupled quantum field theories, the anti-de Sitter/conformal field 
theories (AdS/CFT) correspondence lH3| states that a d-dimensional weakly coupled dual gravitational de- 
scription in the bulk is equivalent to a (d — l)-dimensional strongly coupled conformal field theory on the 
boundary. In recent years, this principle has been used to provide some meaningful theoretical insights in 
order to understand the physics of high T c superconductors from the gravitational dual J, |5j . It was found 
that the spontaneous U(l) symmetry breaking by bulk black holes can be used to construct gravitational 
duals of the transition from normal state to superconducting state in the boundary theory, which exhibits 
the behavior of the superconductor [f|. Due to the potential applications to the condensed matter physics, 
there have been a lot of works studying various gravity models with the property of the so-called holographic 
superconductor (for reviews, see Refs. 7-^j and references therein). 

In most cases, the studies on the holographic superconductors focus on the probe approximation where 
the backreaction of matter fields on the spacetime metric is neglected. When taking the backreaction into 
account, it was found that even the uncharged scalar field can form a condensate in the (2 + l)-dimensional 
holographic superconductor model [5[. Furthermore, in the p-wave holographic dual models, it was argued 
that the phase transition that leads to the formation of vector hair changes from the second order to the first 
order when the gravitational coupling is large enough . There have been accumulated interest to study 
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23j | . Recentl y, t he effect of the backreaction 
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the holographic superconductor away from the probe limit 
has been investigated between holographic insulator and superconductor 

Almost all works on the holographic dual models away from the probe limit were based on numerical 
computations. In order to back up numerical results and gain more insights in the effect of the backreaction, 

29, 3fl, Kanno 



a fully analytic study is called for. Refining the analytic matching method developed in 
calculated the critical temperature of (3 + l)-dimensional holographic superconductors in Einstein-Gauss- 
Bonnet gravity with backreaction and found that the backreaction makes condensation harder [3l|. This 
analytic approach has been extended to derive the critical magnetic field in holographic superconductors with 
backreaction 



32(| . However, the analytic matching method can keep yalid only when the matching point is 
chosen within an appropriate range in higher dimensions {d > 6) [30(. Moreover, when the scalar mass is 
zero in the Gauss-Bonnet holographic superconductors, the curvature correction term does not contribute to 



the analytic approximation, which leads the analytic procedure to break down [30j . Recently, Siopsis and 
Therrien developed a new analytic method. They extended the variational method for the Sturm-Liouville 
(S-L) eigenvalue problem to analytically calculate the critical exponent near the critical temperature and 
found that the analytical results obtained by this way are in good agreement with the numerical findings 
33)]. Considering the effectiveness and accuracy of the S-L method, many authors have us ed it to analytically 
investigate th e prop erties of holographic superconductors in AdS black hole backgrounds 



backgrounds 
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J lj and soliton 



■45j. But these attempts were limited in the probe limit. It was argued that the S-L method 
is more effective for the analytic study of the condensation than the matching method [42]. It is of interest 
to examine whether the S-L method is still valid to explore the holographic superconductivity when the 
backreaction is turned on. This is not trivial since the analytic study can help to confirm the numerical 
result. Furthermore it can clearly disclose the influence of the role of the backreaction in the condensation. 
In this work, we will generalize the variational method for the S-L eigenvalue problem to study holographic 
superconductor away from the probe limit. 

The organization of the work is as follows. In Sec. II, we will introduce the holographic superconductor 
models with backreactions in the d-dimensional AdS black hole background. In Sec. Ill we will give an 
analytical investigation of the holographic superconductors by using the S-L method. We will conclude in the 
last section of our main results. 

II. HOLOGRAPHIC SUPERCONDUCTOR MODELS WITH BACKREACTIONS 

We begin with the general action describing a charged, complex scalar field in the d-dimensional Einstein- 
Maxwell action with negative cosmological constant 

1 ,„ 1 



S = J d d x^ 



^{R~ 2A) - | V" - |W - iqA^\ 2 - m 2 |V| 2 



(1) 



where n 2 = 8irGd is the <i-dimensional gravitational constant, A = — (d— l)(d— 2)/(2L 2 ) is the cosmological 
constant, A and ip represent the gauge field and a scalar field with charge q respectively. Since we are interested 
in including the backreaction, we will take the metric ansatz for the (f-dimensional planar black hole 

dr 2 

ds 2 = -f{r)e-^ r Ut 2 + J^+ r 2 h l3 dx l dx\ (2) 

where / and x are functions of r only, hijdx l dx J denotes the line element of a (d— 2)-dimensional hypersurfacc 
with the curvature k = 0. The Hawking temperature of this black hole, which will be interpreted as the 



temperature of the CFT, is given by 



T H = 



f (r + )e-x< r +)/ 2 
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(3) 



where the prime denotes a derivative with respect to r. r + is the black hole horizon determined by f{r+) = 0. 
We consider the electromagnetic field and the scalar field in the forms 



A = 4>(r)dt, ip = tp{r), 



(4) 



where without loss of generality ip(r) can be taken to be real. Thus, from the variation of the action with 
respect to the matter and metric we obtain the equations of motion 

, 4n 2 r ( , 2 q 2 e x (p 2 ip 2 



X 



r- 



(d-l)r (d-3)/ 
L 2 r 



d-2 



2n 2 r 



P 



0. 



d-2 



d-2 



99 1^/9 f,o q 2 e x 6 z ip z 
roV + -e x <P' 2 + / U' 2 + q J 2 V 



= 0, 



2q 2 ^ 2 



= 0, 



d-2 x ' f\ ™ 2 , , q 2 e*4> 2 



■ip = 0. 



(5) 



(6) 



(7) 



(8) 



It should be noted that the transformation <j) = <fi/q and tp = ipjq'vn. the action ([!) does not change the 
form of the Maxwell and the scalar equations, but the gravitational coupling in the Einstein equation changes 
k 2 — > K 2 /q 2 . Thus, the probe limit is equivalent to letting q — > oo. Without loss of generality, we can set 
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5=1 and keep n 2 finite when we take the backreaction into account 

For the normal phase, tp(r) = 0, we find that x is a constant and the analytic solutions to Eqs. ([6]) and (O 
lead to the AdS Reissner-Nordstrom black holes with the metric coefficient 



f = il J_ 

J T 2 „d-3 
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where \x and p are interpreted as the chemical potential and charge density in the dual field theory respectively. 
When k = 0, the metric coefficient / goes back to the case of the Schwarzschild AdS black hole. 

In order to get the solutions in superconducting phase, where ip(r) =^ 0, we have to count on the appropriate 
boundary conditions. At the horizon r+, the metric functions x an d / satisfy 



X'(r+) 
f(r + ) 



4«; 2 r + 



^'(r + ) 2 + 



d-2 
(d-l)r + 2n 2 r + 



L 2 



d-2 



e^ r +^(p'(r + ) 2 tp(r + f 
m 2 iP{r+) 2 + ^e x ^4>'(r+f 



(10) 



and the regularity condition gives the boundary conditions 
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At the asymptotic AdS boundary (r — > oo), the asymptotic behaviors of the solutions are 



- /J - 



r d— 3 
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(11) 



(12) 



where the exponent A± is defined by [(<i— 1) ± y/ (d — l) 2 + 4m 2 ]/2. Notice that, provided A_ is larger than 
the unitarity bound, both ip_ and tp+ can be normalizable and they can be used to define operators on the 
dual field theory, ip_ =< O- >, tp+ =< C+ >> respectively J, [5j. For simplicity, we will scale L = 1 in the 
following calculation. 



III. ANALYTICAL INVESTIGATION OF THE HOLOGRAPHIC SUPERCONDUCTORS 



Here we will apply the S-L method [33j to analytically investigate the properties of the s-wave holographic 
superconductor phase transition with backreactions. We will derive the relation between the critical temper- 
ature T c and charge density p near the phase transition point and examine the effect of the backreaction. 

Introducing a new variable z = r + /r, we can rewrite the Einstein, Maxwell and the scalar equations into 



X 



Ak 2 
d-2 



Z 3 f 2 



0. 



/ 



, (d-3)/ , (d-l)r 



L 2 z 3 



2n 2 r 2 + 
{d-2)z 3 



1 



= 0, 



X' d-A 
2 z 



2r 2 V 2 

Z f 



d-4 f 
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(13) 
(14) 
(15) 
(16) 



/ J T z A V / P 
where the prime now denotes the derivative with respect to z. 

Since the value of the scalar operator < 0+ > (or < 0_ >) is small near the critical point, we can introduce 

it as an expansion parameter 



€=<Oi >, 



(17) 



with i = + or i = — . Note that we are interested in solutions whe re ib is small, therefore from Eqs. (fT"5|) and 



TBI we can expand the scalar field ip and the gauge field tfi as 
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V> = eipx + e 3 ip 3 + e 5 7/) 5 + 



= 4>a + e 2 02 + e 4 ^4 



(18) 



where e <C 1. The metric function f(z) and x( z ) can be expanded around the Reissner-Nordstrom AdS 
spacetime 



/ = fa + e 2 f 2 + e 4 U + 



X = e 2 X2 + e 4 X4 + 



(19) 



For the chemical potential fx, we will allow it to be corrected order by order 46 1 



(20) 



where <5/i2 > 0. Thus, near the phase transition, we find a result for the order parameter as a function of the 
chemical potential 

1/2 



i« - Ik) 



S^2 



(21) 



whose critical exponent /3 = 1/2 is the universal result from the Ginzburg-Landau mean field theory of phase 
transitions. Obviously, the order parameter becomes zero and phase transition can happen if fi — > /iq, which 
shows that the critical value of /x is ji c = jiQ . 

At the zeroth order, we can get the solution (j>o from Eq. (fT5]) , i.e., the electromagnetic field behaves like 
4>o(z) — /io(l — z d ~ 3 ), which gives a relation /io = p/ r + 3 - At the critical point /z c , we can find [1q — fi c — 
p/r^~ 3 , where r+ c is the radius of the horizon at the critical point. In order to use the analytical S-L method 
33j | . we will set 



,(z) = Ar +c (l-z d - 3 ) 



with A = p/r +c . Inserting this solution into Eq. ([M)) . we obtain the metric function 



(22) 



f (z) = rlg(z)=r 2 



^_z^_ (d-3)^X\ d _ ' 
L 2 z 2 L 2 d-2 [ ' 



(23) 



where we define a new function g[z) for simplicity in the following calculation. 

At the first order, the asymptotic AdS boundary conditions (z — > 0) for ip can be expressed as 



(24) 



So we introduce a trial function F{z) near the boundary z = 



33] 



(25) 



where we have imposed the boundary condition F(0) = 1 and -F'(O) = 0. Substituting Eq. (|2"5l into Eq. p^|) , 
we obtain the equation of motion for F(z) 



F 



2(A, + 2)-rf f/_ 
z 9 



F' 



At {Ai+3-d , g'\ , A 2 (l-z d - 3 ) 2 m 

z 
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z 4 g 2 



z 4 g 



F = 0. 



In order to simplify the following calculation, we will express the backreacting parameter n as 



K n = nAn, n = 0, 1, 2, • 



(26) 



(27) 



where An = — k„ is the step size of our iterative procedure. Considering the fact that k 2 A 2 = k 2 A 2 = 
K n(^ 2 U„-i) + 0[(^ K ) 4 ] (note that we have set k_i = and X 2 \ K _ 1 — 0), we will use the following form of g{z) 
in our discussion 
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( ., . 1 z d -* (rf-3)^(A 2 U„_J d _ 3 _ d _ 3 
lj ~L 2 z 2 L 2 d-2 1 J ' 



where A 2 | tin _ 1 is the value of A 2 for n n -i- After defining a function which obeys 



T(z) = z^+%d - 2){z^ - 1) - (d - 3)L 2 k 2 (A 2 | K „_J(1 - z^)] 



(28) 



(29) 



we can convert Eq. (J26J) to be 
(TF'Y + T 



A t (Ai + 'i-d , g'^ , A 2 (l-2 d - 3 ) 2 m 

z 
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z 4 9 2 



z 4 g 



F = 0. 



(30) 



From the Sturm-Liouville eigenvalue problem [47| , we write down the expression which can be used to estimate 
the minimum eigenvalue of A 2 

-_fiT{F»-UI*)dz 



IoTVF*dz 



(31) 



with 



u= Aj, / \ + 3-d | g' 

z \ z g 



z 4 g 



V = 



(1 - z d ~ 3 ) 2 
zW ■ 



(32) 



In order to use the variation method, we will assume the trial function to be F{z) = 1 — az 2 , where a is a 
constant. 

Using Eq. (1311) to compute the minimum eigenvalue of A 2 for i = + or i = — , we can obtain the critical 
temperature T c for different strength of the backreaction k and the mass of the scalar field m from the following 
relation 



Tr 



1 

-in 



X, 



(33) 



As an example, we calculate the case for d = 5 and m 2 L 2 — —3 with the chosen values of the backreaction 
parameter k for i — +, i.e., A + = 3. Setting An = 0.05, for kq = we have 



2(-18 + 27a- 14a 2 ) 
6(4 In 2 - 3) + 16(3 In 2 - 2)a + (24 In 2 - 17)a 2 



^ = n ^T^^TT^n o 2 . (34) 



whose minimum is A 2 | Ko = 18.23 at a = 0.7218. According to Eq. ([33| . we can easily get the critical 



temperature T c = 0.1962/j 1 / 3 , which agrees well with the numerical result T c = 0.1980/9 1 / 3 
substituting A 2 | reo into Eqs. (j2"8"j) and (|2T))) we obtain 



. For ki = 0.05, 



A 2 _ 4.466- 6.682a + 3.462a 2 

~~ 0.1714 - 0.1600a + 0.04592a 2 ' ' 

which attains its minimum A 2 | Kl = 18.11 at a = 0.7195. Hence the critical temperature reads T c = 0.1934/? 1 / 3 , 

which is also in good agreement with the numerical result T c = 0.1953/5 1 / 3 . For K2 = 0.10, putting A 2 | Kl in 

Eqs. (HHJ) and (JUJ) we arrive at 

A 2 _ 4.364- 6.478a + 3.349a 2 

~ 0.1740 - 0.1633a + 0.04705a 2 ' ( ' 

whose minimum is A 2 | K2 = 17.75 at a = 0.7122. So the critical temperature is T c = 0.1852/9 1 / 3 , which is again 

consistent with the numerical finding T c = 0.1874/5 1 / 3 . For other values of k, the similar iterative procedure 

also can be applied to present the analytic result for the critical temperature. When we reduce the step size, 

for example to fix Ak = 0.025, we can also compute the critical temperature T c in the similar way. In Table 

U we give the critical temperature T c for the scalar operator < C + > when we fix the mass of the scalar 

field m 2 L 2 = —3 for different strength of the backreaction by choosing the step size Ak = 0.05 and 0.025, 

respectively. We find that the analytic results derived from the S-L method are in very good agreement with 

the numerical calculation. Furthermore we observe that when we reduce the step size Ak, we can improve 

the analytic result and get the critical temperature more consistent with the numerical result. 

TABLE I: The critical temperature T c with the chosen values of the backreaction parameter k and the step size Ak 
for the condensates of the scalar operator < O+ > in the case of 5-dimensional AdS black hole background. Here we 
fix the mass of the scalar field by m 2 L 2 — —3. 

Analytical(AK = 0.05) Analytical(AK = 0.025) Numerical 

k = 0.1962p 1/3 0.1962p 1/3 0.1980p 1/3 

k = 0.05 0.1934p 1/3 0.1934p 1/3 0.1953p 1/3 

k = 0.10 0.1852p 1/3 0.1853p 1/3 0.1874p 1/3 

k = 0.15 0.1718p 1/3 0.1722p 1/3 0.1748p 1/3 

k = 0.20 0.1540p 1/3 0.1549p 1/3 0.1580p 1/3 

k = 0.25 0.1330p 1/3 0.1345p 1/3 0.1382p 1/3 

k = 0.30 0.1098p 1/3 0.1123p 1/3 0.1165p 1/3 



For completeness, we also extend the investigation to the 4-dimensional AdS black hole background. In 
Table UH we present the critical temperature T c of the chosen parameter k with the scalar operators < 0_ > 
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and < 0+ > for the (2 + l)-dimensional superconductor if we fix the mass of the scalar field by m 2 L 2 = — 2 
and the step size by Ak — 0.05. The agreement of the analytic results derived from S-L method with the 
numerical calculation shown in Tables HI and IPTl is impressive. 



TABLE II: The critical temperature T c obtained by the analytical S-L method (left column) and from numerical 
calculation (right column) with the chosen values of the backreaction parameter k for the condensates of the scalar 
operators < O- > and < 0+ > in the case of 4-dimensional AdS black hole background. Here we fix the mass of the 
scalar field by m 2 L 2 — —2 and the step size by Ak = 0.05. 







< O- 


> 


<o+ 


> 


K 


= 


0.2250p 1/2 


0.2255p 1/2 


0.1170p 1/2 


0.1184/9 1 / 2 


K 


= 0.05 


0.2249p 1/2 


0.2253p 1/2 


0.1163p 1/2 


0.1177/9 1/2 


K 


= 0.10 


0.2246p 1/2 


0.2250p 1/2 


0.1141p 1/2 


0.1156/9 1/2 


K 


= 0.15 


0.2241p 1/2 


0.2245p 1/2 


0.1106p 1/2 


0.1121/9 1/2 


hi 


= 0.20 


0.2235p 1/2 


0.2239p 1/2 


0.1057p 1/2 


0.1074 / 9 1/2 


hi 


= 0.25 


0.2226p 1/2 


0.2230p 1/2 


0.0998p 1/2 


0.1017p 1/2 


K 


= 0.30 


0.2216p 1/2 


0.2220p 1/2 


0.0929p 1/2 


0.0951/9 1/2 



It further supports the observation obtained first in the nu mer i cal 
action can make the scalar hair more difficult to be developed 
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computation that the stronger backre- 



23| . The consistency between the analytic 



and numerical results indicates that the S-L method is a powerful analytic way to investigate the holographic 
superconductor even when we take the backreaction into account. 

IV. CONCLUSIONS 

We have generalized the variational method for the Sturm-Liouville eigenvalue problem to analytically 
investigate the properties of the holographic superconductor with backreactions. We found that in the fully 
backreacted spacetime, the S-L method is still powerful to disclose the property of the condensation. Our 
analytic results are in very good agreement with those obtained from numerical computations. If we reduce 
the step in the iterative procedure, we can further improve our analytic results and improve the consistency 
with the numerical findings. Our analytic result shows that the backreaction makes the critical temperature 
of the superconductor decrease, which can be used to back up the numerical finding as shown in figure 2 of 
Ref. 5] that the backreaction can hinder the condensation to be formed. 
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